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Low-energy properties of the one-dimensional Kondo lattice model are investigated by 
using bosonization and renormalization group methods. The formation of the spin and 
charge excitation gaps at half-filling is discussed in detail both for antiferromagnetic and 



ferromagnetic Kondo couplings. Also, the effect of the interaction between conduction 
electrons on the gap formation is studied both for the Kondo insulating phase and 
for the Haldane gap phase. It is found that while the repulsive interaction changes 
the dependence of the excitation gap on the Kondo exchange coupling, the attractive 

interaction can drive the Kondo (or Haldane) insulating phase to a completely different 

<N ■ 

^^ ■ phase of massless metallic states. 
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q §1. Introduction 



The Kondo lattice model is a fundamental model for heavy electron systems such as Cerium 
or Uranium compounds. One-dimensional (ID) version of the model has been studied ex- 
tensively for its simplicity by using numerical as well as analytic methods.a'SSoSBQ'BtP 
The model shows various phases according to the electron density and the strength of the 
Kondo coupling. At half-filling the Kondo insulating phase realizes, which has the excita- 
tion gaps both in the spin and charge sectors.H'EP Away from half-filling the model shows the 
Tomonaga-Luttinger liquid phase for small Kondo exchange couplings.Q'&cP More recently, 
the effect of interaction between conduction electrons on these phases has been studied nu- 
merically.li3 ) 



We have previously investigated this model by using bosonization methoao, which is a 
powerful tool to study low-energy properties of ID quantum systems. In this method, one 
needs to linearize the dispersion of electrons at the first stage for deriving a low-energy ef- 
fective hamiltonian. However, since localized spins have no dispersion, it is rather difficult 
to apply this method to the Kondo lattice model straightforwardly In the previous paper, 
to bypass this difficulty we naively introduced a finite velocity for the excitation of localized 
spins in the microscopic model, and then performed bosonization.© This prescription, how- 
ever, may not give desirable functional forms of the spin and charge gaps for the original 
Kondo lattice model: the predicted gaps are not consistent with those obtained by numer- 
ical studies.® In this paper, by starting with the original Kondo lattice model, we develop 
a low-energy effective theory improved in this respect, and reexamine the formation of the 
spin and charge gaps at half filling both for antiferromagnetic (AFM) and ferromagnetic 
(FM) Kondo couplings. Based upon this formulation, we also study how the short-range 
interaction between conduction electrons affects the spin/charge gap formation. 

The organization of the paper is as follows. In §2 we introduce the bosonized effective 
hamiltonian for the ID Kondo lattice model, and derive the scaling equations for various 
couplings which characterize the low-energy fixed point properties. In §3 we first discuss 
the case of the AFM Kondo coupling, and describe the formation of the spin and charge 
gaps in the Kondo insulating phase. The effect of interaction between conduction electrons 
on this phase is also investigated. It is shown that while the repulsive interaction destroys 
neither the spin nor charge gaps, the strong attractive interaction should change the Kondo 
insulating phase to a new gapless phase characterized by the Luther-Emery model and the 
s = 1/2 spin chain model. In §4 the case of the FM Kondo coupling is also discussed in 
detail. Conclusion is given in §5. 

§2. Bosonization and scaling equations 

We first introduce the model hamiltonian for the ID Kondo lattice, and apply the bosoniza- 
tion method to obtain the scaling equations which describe the fixed point properties. We 
consider the following hamiltonian, 

H = -£^ct,cwi + h.c. + Uj2 c \i c ^ iC \i c ^ + J kY1 c li2 (T <xP c fc ' S fi> ( 2 ' X ) 

i,cr i i,aj3 

where c a ; L (c ai ) is an annihilation (creation) operator for conduction electrons at z-th site, 
and Sfi is an operator for localized s — 1/2 spins. The second term is the on-site interaction 
between conduction electrons, and the last term is the Kondo exchange interaction. Note 
that we have not introduced the RKKY interaction here, which was added to the microscopic 
hamiltonian in our previous paper & In order to apply renormalization group methods to this 



model (2T), we first bosonize the charge and spin degrees of freedom of conduction electrons 
and take their continuum limit with preserving full degrees of freedom of localized spins at 
each lattice site. Applying abelianli^ and non-abelianliS t3* bosonizations to charge and spin 
degrees of freedom of conduction electrons respectively, we have the effective hamiltonian, 

H = H C + H c . f} (2.2) 

H -= /dx[-^(40c(x)) 2 + ^n2(x)] + ^ fdxcos(V^(f> c (x) + 5x) 
J zA c I a. J 

+ / dxv s [J cL (x) ■ J cL (x) + J cR (x) ■ J c r(x)} - X dxJ cL (x) ■ J cR (x), (2.3) 
H c _ f = J KA I dx[J cL (x) + J cR {x)] ■ S f (x) 

+Jkb [ — [e 2ikFX -tr(g(x)(r)e lV ^ Mx) + h.c] ■ S f (x), (2.4) 

J a 2 

where 4> c (x), n c (x) are boson fields for the charge degrees of freedom of conduction electrons 
and its conjugate momentum, and J c l,r is a spin current operator which satisfies level-1 
SU(2) Kac-Moody algebra, and g(x) = c aL cp R (a, (3 =|, j) is a fundamental representation of 
SU(2). Sf(x) = J2i fai( a af3/^)f(3iS(x — Xi) is the spin density operator for localized spin. The 
U um - and A-terms are, respectively, the Umklapp and backward scattering terms of electron- 
electron interaction, and the parameter 5 defined by 5 = Akp— 2ir measures the deviation from 
the half-filling case. The forward scattering between conduction electrons is incorporated in 
the Tomonaga-Luttinger parameter K c , via the initial value K c = 1/ wl + 2U/ / Kv-p. 

In the microscopic model ( |2.1| ), the bare parameters should read Jka = Jkb = Jk an d 
U um = X = U, which are to be modified via the renormalization process. Note that if 
we regard the initial values of K c , U um , and A as independent parameters, this effective 
hamiltonian describes the low-energy properties of a more general model than eq. (|2.1| ), which 
has momentum-dependent interactions. We will also extend our arguments to such a general 
model in the subsequent sections. In what follows we mainly consider the half-filling case 
5 = 0. Thus we impose the commensurate condition, 2kpa = n where a is the lattice 
constant. 

Now we derive the scaling equations for various coupling constants of the effective hamil- 
tonian. For this purpose we utilize current algebra techniques for the operators. The charge 
field, C , satisfies the operator product expansion (OPE) of U(l) Gaussian model,E-$ 

_io^(*,T)_-ia« c (0,0) ta AMW) , 9- z (j) c (0,0) 
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with z = x + iv c r, z = x — iv c r. J c l,r and g(x) satisfy the OPE of level-1 SU(2) Wess- 
Zumino-Witten modelPBS 

r p-abc 

•&(*, r)J b cL (0, 0) ~ 4 + — J C C L (0, 0) + • ■ -, (2.7) 

z z z 

J c a L (x,r)^(0,0)~^(0,0) + ---, (2.8) 

tr[g(x,T)a a }tr[g\0,0)a b ] ~ ^ + Q V bc J c c L (0, 0) + g) V 6 V c c fl (0, 0) 

\A/\ As A> 

+|y[^(0, 0) J c a L (0, 0) + ^ 2 J^(0, 0) J c yo, 0)] + • • , (2.9) 

with z = x + w> s r, = x — if s r, and t a , the generator of the SU(2) Lie algebra. We omit 
those terms which are irrelevant to the following arguments. Using these formulae, we can 
derive scaling equations, 

dJKA = JjA + JjB h{ ^ (21Q) 

alnL irv s irv s v s 
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In the dilute limit of localized spins, the velocity of conduction electrons is not renormalized. 
In this case, if one neglects backward and Umklapp processes of electron-electron interaction, 
and expand the Luttinger liquid parameter up to the lowest order in the forward scattering, 
K c ~ 1 — U/tt, one can reproduce the scaling equations for the single impurity Kondo problem 
in Luttinger liquids, which were obtained by g-ology methods.^ In the lattice system, the 
RKKY interaction between localized spins is generated by higher order processes of the 
Kondo exchange interaction. With the use of eq.( |2.9| ), we obtain the RKKY interaction 
from the second order of Jkb, 

H RKKY = -J 2 KB fdx fdx'f 3 (^) COS2kF [ X ; X,) S f (x) ■ S f (x'), (2.17) 

<J <J (J Q iXj <Xj 



with 
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/a(«) = / de ~— Q ; ■ 2m ^ - ( 2 - 18 ) 

Jo (cos 2 6» + w 2 sm 2 6»)~ 

The predominant contribution comes from the nearest neighbor interaction. Since 2kpa = it, 

the RKKY interaction is antiferromagnetic for nearest neighbor sites. Thus we replace the 

RKKY interaction term approximately by 

Hrkky = Jkb^LJ^—) ^2 Sfi ' Sfi+i- (2-19) 

v s i 

Then the scaling equation for the RKKY interaction is given by 

dJKKKY = Jk /8( «E). (2 . 20 ) 

a in L ttv s v s 

Eqs. fl2.10|) ~( |2.l4D and (|2.20 ) determine the low-energy properties of the model. 



We note here that in our previous paper the RKKY interaction was naively introduced in 
the microscopic hamiltonian.cP In the following we will take into account that the RKKY 
term is a higher-order term generated by the Kondo interaction, and see what follows from 
this fact. 

§3. Antiferromagnetic Kondo coupling case 

We start our discussions with the case of AFM Kondo coupling (Jk > 0). In this case, 
from eqs.( p.lOJ ) and (|2.20Q , we see that Jka and Jrkky how to the strong-coupling regime. 



Since initially Jk > Jrkky — 0, the Kondo interaction Jk is predominant over the RKKY 
interaction Jrkky- Thus low-energy properties at the fixed point are determined essentially 
by Jk- We have three different situations according to the interaction strength for conduction 
electrons: case A: U = and K c (0) = 1, case B: U > and K c (0) < 1, and case C: U < 
and K c (0) > 1. 

3.1 case A: U = and K c (0) = 1 

In the absence of interaction between conduction electrons, U = 0, the scaling equations 
for the Kondo couplings are simplified as, 

dJ KA __J 2 KA + JlB h{ v_ 1 ^ (31) 



dhiL tiv s irv s " v. q 
dJxB 2JkaJkb 



(3.2) 
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Thus Jka and Jkb are marginally relevant, and produce spin and charge gaps. The fixed 
point is the Kondo insulator as pointed out beforeBB 1 If one neglects the renormalization of 
the velocity, the spin gap A s and the charge gap A c generated by Jka and Jkb are given by 

A SiC ~e"^, (3.3) 



for J K C?; s . The Jk- dependence of the spin gap coincides with that obtained by Tsunetsugu 
et all. 

3.2 case B: U > and K c (0) < 1 

When the on-site interaction for conduction electrons is repulsive, the resulting A-term 
(backward scattering) is marginally irrelevant, and negligible. Jkb is a relevant interaction 
with scaling dimension 3/2 + K c /2 as seen from eq. (|2.11|) . Since this dimension is larger 
than that of the Umklapp interaction, the Umklapp interaction U um grows faster than Jkb 
at half-filling, which may give rise to the Mott-Hubbard gap in the charge sector of conduction 
electrons. However if the initial value of U um is much smaller than that of Jkb, Ji^s-term 
may generate the charge gap before U um -term flows into the strong-coupling regime. Thus 
several different cross-over behaviors are expected to appear according to the initial values 
of parameters. In the following we show the expressions for spin and charge gaps in three 
limiting cases. 

« '-^ » 2 ~^W and U U0)/Jk(0)>1- 
The Mott-Hubbard gap due to C/ um -term opens prior to the spin gap formation due to J K 



in this case. After the charge gap formation, K c = 0, and eq. (|2.11| ) becomes, 

dJKB Jkb 



(3.4) 



dlnL 2 
Thus the spin gap for small Jk is given by 

A s ~ J K - (3.5) 

The charge gap is given by the Mott-Hubbard gap of the ID Hubbard model, if Jk is 
sufficiently small compared to U. 

(ii) i^M » ^£0 and Uum (0)/J K (0) « 1. 

In this case, the charge gap generation due to J^s-term is predominant over the Mott 
transition due to the Umklapp interaction between conduction electrons. Thus the spin and 
charge gap is given by 

2 



&s,c~Jk KC : (3-6) 

for small Jk- The realization of this case requires a large forward scattering of electron- 
electron interaction and a small Umklapp interaction at the initial stage of the renormal- 
ization process. However, in the microscopic model given by eq. (|2.1| ), these parameters are 
initially equal to each other. Thus in order to actually observe this behavior, we should con- 
sider a more general model such as the g-ology model which has different coupling constants 
of electron-electron interaction depending on transfered momentum.liS' 



(■■■\ 1-K40) „ 2J KA (0) 

V 11 ^ 2 ^ TTVs(O) 



In this case, the first term of the right hand side of eq.( |2.11 ) is negligible. Thus the spin 



gap has the same J^-dependence as the case of K c = 1 (U = case). 

TTV S 

A s ~ e j k . (3.7) 

The charge gap in this case is generated by U um -teini or J^s-term depending on the initial 
value of these couplings. 

It is to be noted here that the above different behaviors observed in the gap formation 
do not imply that there appear new different universality classes: the system still belongs 
to the same universality of the Kondo insulating phase as far as U > 0, and cross-over 
behaviors rather than phase boundaries are seen when the initial values of bare parameters 
are changed. 

We would like to compare the above results with those obtained by Shibata et al.t^P. 
They calculated numerically the spin gap of the model given by eq.( p.l|) . They concluded 
that the spin gap is given by A s ~ exp(— TXV a /aJjc) for small Jk and any U, where a 
parameter a depends on U. Their conclusion does not seem to coincide with our results 
in the case (i) shown above. Actually, the parameter region where their calculation was 
performed corresponds to the case (iii) or the case where (1 — K c (0))/2 is of the same order 
as Jka(0)/ttv s (0). We have solved numerically the scaling equations when (1 — K c (0))/2 is 
of the same order as Jka(0)/ttv s (0), and deduced the functional dependence of the spin gap 
on Jk- We found that the results fit in with the exponential function fairly well. Thus our 
results are consistent with those of Shibata et al. We think that the value of Jk they used 
may not be small enough to observe the power-law dependence of eq.( 



3.3 case C: U < and K c (0) > 1 

In the attractive case of U, the Umklapp term is irrelevant, and can be discarded. However 
since A < 0, the A-term is marginally relevant. The Jxs-term may be relevant or irrelevant 
depending on K c (0) and Jk- Therefore we have three different behaviors as shown below. 

(i) J kb predominates over Jka and A. 

If \U\ is sufficiently small compared to Jk, and K c (0) is not so large, the J^s-term is still 
relevant. The spin and charge gaps are generated by the Jka- and J^s-terms. For small 
Jk, the gaps are given by 

A s , c ~e"^. (3.8) 

Note that the universality class still belongs to that of the Kondo insulator, 
(ii) Jka predominates over Jkb and A. 



If K c (0) is sufficiently large, the J^s-term becomes irrelevant. Moreover if |A(0)| is smaller 
than J K , J K a grows to a strong-coupling value before A flows into the strong-coupling regime. 
In this case, the J^A-term generates the spin gap given by 

A s ~e"^?. (3.9) 

The charge sector is massless. Thus the fixed point of the charge degrees of freedom is the 
Tomonaga-Luttinger liquids described by c = 1 U(l) Gaussian model. This would bring 
about a new phase driven by the interaction. However, we recall that for the microscopic 
model given by eq.( |2.1| ), we can not control the parameters K c (0) and A(0) independently. 
Thus it is not obvious whether this phase may realize for the model ( |2.1| ). If one considers 
the model in which the strength of the forward scattering of electron-electron interaction 
which controls K c (0) and the backward scattering A(0) are independent parameters, one can 
surely find this phase. 

(iii) A predominates over Jka and Jkb- 

We now discuss the case where \U\ is sufficiently large and \U\ ^> Jr- Then the Jkb- 
term is irrelevant, and thus the charge sector becomes massless, forming the Tomonaga- 
Luttinger liquid described by c = 1 U(l) Gaussian model. On the other hand, the spin gap 
generation occurs due to the A-term (backward scattering term) prior to that due to the 
Jft-^-term. This state obtained for conduction electrons is nothing but a metallic state with 
the spin gap which is often referred to as the Luther-Emery class. Localized spins are now 
completely decoupled from conduction electrons, making a massless mode classified as SU(2) 
Tomonaga-Luttinger liquid described by c = 1 level- 1 SU(2) Wess-Zumino-Witten model. 
Therefore it turns out that a sufficiently large attractive interaction should drive the Kondo 
insulator to a completely different phase, in which both of the spin and charge sectors are 
massless: conduction electrons are in the Luther-Emery phase and localized spins are in 
SU(2) Tomonaga-Luttinger liquid. Recently Shibata et al. have also found this massless 
phase for large \U\ by using the density matrix renormalizatoin group method.B°P 

§4. Ferromagnetic Kondo coupling case 

We now move to the discussions for the case of FM Kondo couplings. It is known that the 
spin gap a la Haldane opens at half filling in this caseBElP We now explore how the formation 
of the spin and charge gaps is affected by the electron-electron interaction for conduction 
electrons. 

4.1 U = andK e (Q) = 1 
In this case, as is seen from eqs. (|2.10|) and (|2.11| ), the Jka- and J^s-terms are marginally 



irrelevant, and these couplings scale to zero. In the process of the renormalization, Jrkky 



becomes larger than \Jka\ and \Jrb\- Thus we can not discard the RKKY interaction 
compared to the Kondo coupling. The induced RKKY interaction gives rise to a dispersion 
of localized spins, and thus changes the scaling equations obtained in §2. This point is quite 
contrasted to the antiferromagnetic case discussed above, where the Jra- and J^-terms 
are marginally relevant. For Jrkky 3> |^a|, \Jkb\, we can bosonize localized spins in the 
continuum limit, and treat the Kondo coupling perturbatively. Then the field theoretical 
limit of the RKKY interaction is 

Hrkky = J dxv f [J fL {x) ■ J fL {x) + Jfn{x) ■ J jr{x)] - X f J dxJ fL {x) ■ J fR {x), (4.1) 

where J/l,r satisfy level-1 SU(2) Kac-Moody algebra. Xf > is a coupling constant for 
the marginally irrelevant operator which exists because of chiral SU(2)xSU(2) symmetry. 
This is level-1 SU(2) Wess-Zumino-Witten model with marginally irrelevant interaction. We 
now incorporate this RKKY term into a non-perturbed part of the hamiltonian. Then the 
bosonized expression of the Kondo exchange interaction is altered to 

H c - f = JraJ dx(J cL (x) + J cR {x)) ■ (J/l(x) + J fR {x)) 

+J KB J dx[tr(g(x)^) ■ tr(g f (x)^)e^^ + h.c], (4.2) 

where gj(x) = f a Li x )fpR{ x ) ( a i/3 =T>!)- The scaling dimension of the second term of 
eq. (|4.2j ) is now 1 + K c /2. Note that K c is renormalized by the J^s-term and should be 
smaller than 1. This term is relevant and generates the spin and charge gaps. The low-energy 
fixed point is nothing but the Haldane gap phase.cf Although the spin and charge gaps are 
given by the power-law function of Jkb, we cannot write down their expressions explicitly 
in terms of bare coupling J^, because Jkb i n eq.( [4.2| ) is renormalized, being different from 
bare Jk- 

4.2 U > and K C (Q) < 1 

If we introduce the repulsive interaction for conduction electrons, Jkb becomes relevant, 
as is seen from eq. fl2.lip . Thus the spin and charge gaps open as discussed in §3.2. The spin 
gap is given by eq. (|3.5| ) or ( p.6|) . The universality class is still characterized by the insulating 
phase with the Haldane- gap. 

4.3 U < andK c (0) > 1 

In the case of attractive interaction, the Jka~ and J^s-terms of eq. fl2.4j) and the Umklapp 
interaction U um are irrelevant. Thus we should take into account the RKKY interaction 
again as in the case of U = 0. We bosonize localized spin and use eq. (|4.2| ) as the effective 



Kondo exchange interaction. Then the scaling equation for Jkb is changed to 



amL v 2 irv s itv s irvf irv c ' 

The J^^-term may be relevant or irrelevant depending on the initial values of A = Ao and 
K c = K c0 , which are renormalized to the values different from U and K c (0). Moreover, 
as is discussed in §3.3, the A-term flows to a strong-coupling regime. Thus we have two 
possibilities: (i) Jkb predominates over A, (ii) A predominates over Jkb- 

(i) Jkb predominates over A. 

If | Ao | is sufficiently small, and K c0 is not so large, Jkb may be relevant. Since Jkb grows 
faster than A, the J^s-term generates spin and charge gaps prior to the spin gap formation 
due to the A-term provided that the initial value of Jkb is not much smaller than Ao- The 
fixed point is the same class as the Haldane gap phase. 

(ii) A predominates over Jkb- 

For sufficiently large negative U, A flows to a strong-coupling regime, A — > — oo, faster than 
Jkb- Then, the right hand side of eq.( |4.3j ) becomes positive. The J^s-term is irrelevant in 
this case. Thus the charge sector of conduction electrons is massless, and described by c = 1 
U(l) Gaussian model, i.e. the Tomonaga-Luttinger liquids, and the spin gap of conduction 
electrons opens due to the A-term, as has been discussed for the AFM case. Therefore, even 
for the FM case, a sufficiently strong attractive interaction brings the system to a massless 
phase with conduction electrons in the Luther-Emery phase and localized spins in the SU(2) 
Tomonaga-Luttinger phase. 

§5. Conclusion 

We have studied the ID Kondo lattice model by means of bosonization techniques, and 
discussed how the mass gap is generated at half filling. It has been confirmed that the spin 
and charge gaps are indeed generated both for the AFM and FM Kondo couplings, realizing 
the Kondo insulating phase and the Haldane-gap insulating phase. By taking into account 
that the RKKY interaction is a higher order term induced by the Kondo coupling, we have 
shown that the generation of the spin gap occurs in different manners for the AFM and FM 
cases. This is one of the central points to be claimed in the present paper, which was not 
considered in our previous paper .tf The results thus obtained for the spin gap are consistent 
with the previous numerical studies, and may explain why the functional dependence of the 
spin gap on Jk can differ between the AFM and FM couplings.B'' 

We have also investigated how the short-range interaction between conduction electrons 
affects on the formation of these excitation gaps. It has been shown that although the 
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repulsive interaction changes the functional dependence of spin and charge gaps on the 
Kondo coupling, the universality class is still characterized by that for U — 0. On the 
other hand, a sufficiently strong attractive interaction can destroy the spin and charge gaps, 
resulting in a completely different class of massless liquids consisting of the Luther-Emery 
state and the SU(2) Tomonaga-Luttinger liquid state. 

We have seen that the spin gap formation can be described rather well in the present 
approach. As for the charge gap, however, we have not been able to completely resolve the 
discrepancy between the present results and the numerical results, although the qualitative 
agreement was already established concerning whether the mass is generated or not. For 
example, the exponential dependence of the charge gap on Jk obtained for the AFM case 
with U = does not seem to fit in with the J^--linear dependence deduced numericallya- 1 To 
resolve this problem is one of the main issues to be addressed in the future work. 

Finally we wish to note that the present results are applicable not only to the model 
with the Hubbard-like on-site interaction but also to more general models with momentum- 
dependent electron-electron interaction. For such general models, one can find the metallic 
phase with spin gap predicted in §3 for the AFM Kondo model with attractive electron- 
electron interaction. 
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